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ABSTRACT 

LSS 2018, the central star of the planetary nebula DS1, has been 
found to be a double-lined spectroscopic binary with a period of 
8.571 hours. Light variations with the same period have been 
observed in U_, B_, and Vj in the wavelength regions defined by the two 
IUE cameras; and in the strength of the Cl 1 1 4647 emission line. The 
light variations can be accurately predicted by a simple 'reflection' 
effect, and an analysis of the light curves yields the angular 
diameter and effective temperature of the primary, the radii of the 
two stars in terms of their separation, and the inclination of the 
system. Analysis of the radial velocities then yields the masses of 
the two stars, their separation, the distance of the system, the 
absolute magnitude of the primary, and the size of the nebula. 


- 3 - 


I. INTRODUCTION 

Bond (1976) and Grauer and Bond (1983) have initiated a search 
for central stars which show periodic light variations. To date, 
five new objects have been found, four of which have turned out to be 
close binaries showing a very large reflection effect. These are the 
central stars of Abell 41 (P = 0^1 1 32269 ; Grauer and Bond 1983), 

Abell 63 (P = 0^46506918; Bond, Li Her, and Mannery 1978), Abell 46 
(P = 0^4717294; Kurochkin 1980); and Kl-2 (P = 0^6707, Bond 1979). 

The central stars of Abell 46 (V477 Lyr) and Abell 63 (UU Sge) are 
also eclipsing binaries, and the light curve of UU Sge has been 
analyzed by Bond, Liller, and Mannery (1978). The spectrum of the 
central star of Abell 41 has been discussed by Green, Liebert, and 
Wesemael (1984). Because of the lack of fundamental mass and radius 
determinations for central stars, it would be of great importance to 
establish accurate radial velocity curves for these systems and to 
obtain time-resolved observations in the ultraviolet, where most of 
the observable radiation is being emitted. Because of the faintness 
of these objects (B > 15), and their short periods, this has not been 
done. 

Drilling (1983) has discovered 12 new sdO stars, 3 of which have 
subsequently been found to possess planetary nebulae of low surface 
brightness. Because these objects are all brighter than B = 13 and 
had not been observed previously, the writer monitored them for light 
variation during the nights of 1984 April 3 through 10 using the 0.9- 
meter telescope at Cerro Tololo. Only one, LSS 2018, was found to 
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vary, and subsequent time-resolved spectroscopic observations were 
made at high resolution with the 4-meter echelle spectrograph at 
Cerro Tololo and at low resolution with the IUE satellite. These 
observations are described and analyzed for the relevant stellar 
masses and radii below. In all cases, the observations have been 
phased using a period of 8.571 hours with phase zero (maximum light) 
occurring at (heliocentric) J. D. 2445796.671. 

II. OBSERVATIONS 
a) UBV Photometry 

Photoelectric UBV measurements of LSS 2018 were made with the 
0. 9-meter telescope at Cerro Tololo on the nights of 1984 April 3 
through 10 using a standard photometer and filters. The reduction 
procedure outlined by Schulte and Crawford (1961) was used with the 
following mean extinction coefficients: k_ = 0.16, = 0.09, Jc^ = 

-0.03, = 0.31, and = 0.01. Transformation coefficients were 

determined from 3/ observations of standard stars for which 
magnitudes and colors are given by Graham (1982). Night corrections 
determined from frequent observations of a comparison star located 
6^1 E and K6 N of LSS 2018 were applied to all observations of LSS 
2018. The average V^ magnitude and colors of the comparison star are 
V_ = 10.87, B-V_= 0.56, and U-B_ = 0.07. The results for LSS 2018 are 
shown as squares in Fig. 1. Also shown in Fig. 1 (as crosses) are 
some very accurate observations of LSS 2018 made using the B_ filter 
alone. Each point represents the average of as many as 400 five- 
second integrations made during the nights of 1984 April 5 and April 
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10. The scatter of the individual observations about the mean light 
curve is random with a standard deviation of +0.02. 

b) IUE Spectrophotometry 

Low-resolution (6 A) IUE observations were made with both the 
LWP and SWP cameras during 1984 May 20 and 21. These observations 
were reduced by the staff at Goddard Space Flight Center using the 
software described by Turnrose and Thompson (1984). The quantities 
plotted in Fig. 1 are the total fluxes observed in the wavelength 
intervals 1240-1945 A (SWP) and 1945-3120 A (LWP), corrected for a 
reddening of E_(B.-V.) = 0.15 (Schonberner and Drilling 1984) according 
to the reddening law of Seaton (1979). Also plotted as magnitudes 
in Fig. 1 are the FES (fine error sensor) readings from IUE, shifted 
vertically so as to give the best fit to the magnitudes obtained at 
Cerro Tololo. Because the sensitivity function of the FES is inter- 
mediate to those of the Johnson B_ and V_ filters (Holm and Rice 1981), 
and because the variation in B^ over one cycle is negligible, the 
coincidence of the two sets of observations shows that the light 
variations remained strictly periodic during the six weeks which 
elapsed between the photometric observations made at Cerro Tololo and 
those made with IUE. 

c) High-Resolution Spectroscopy 

High-resolution (0.3 A) observations of LSS 2018 giving 70% 
coverage of the wavelength interval from 4600 to 4900 A were made 
with the 4-meter echelle spectrograph and SIT-vidicon detector at 
Cerro Tololo during the nights of 1984 April 13 and 14. The exposure 
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time was 30 minutes in all cases. After removal of the fixed pattern 
noise, the background-subtracted intensities for LSS 2018 were 
divided by those for a quartz comparison lamp, and then corrected for 
the mean atmospheric extinction derived above. The spectra were 
wavelennlh calibrated using ThA comparison spectra obtained 
immediately before and after each observation. 

The flux in the C 1 1 1 4647 emission line, in units of the 
continuum flux in a 3.4 A band underlying this line at minimum light, 
is plotted in Fig. 1. The plotted points (squares) show that the 
strength of the C 1 1 1 4647 emission varies in phase with that of the 
continuum, and that it disappears completely at minimum light. The 
heliocentric radial velocity given by the C 1 1 1 4647 emission line is 
plotted in Fig. 2 (squares), as is that given by the He II 4686 
absorption line (triangles). It is seen that these observations can 
be fit quite well by sine curves of amplitude 148 km/sec and 74 
km/sec, respectively, with phases u and 0, respectively, at maximum 
light. The best fit is obtained for a system velocity of -25 km/sec. 

III. ANALYSIS 

The observations are well explained using the model applied by 
Grauer and Bond (1983) to the central star of Abell 41, in which a 
cool star (from which we see only 'reflected' radiation) orbits a 
very hot subdwarf (which heats both the planetary nebula and the 
facing hemisphere of the cooler companion). The symmetry of the 
radial velocity and light curves are well understood if the orbits 
are circular, the C 1 1 1 4647 emission is produced in the atmospher° of 
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the secondary, and the Hell 4686 absorption is produced in the 
atmosphere of the primary. A small contribution to the Hell 4686 
line by the secondary may be present near maximum light, when the 
contribution of the secondary to the spectrum as a whole is greatest 
accordin'- to this model. 

C 1 1 1 4650-51 (two blended lines) and NIII 4641-42 (near the end 
of one of the echelle orders) are also present in tmission, and both 
the strengths and radial velocities vary in phase with the Cl 1 1 4647 
emission. NV 4603 and 4619 are weakly present in absorption and show 
the same radial velocity variations as Hell 4686. The Hp profile is 
very complex, but appears to consist of the following: (a) an 

absorption component of constant strength which shows the same 
variations in radial velocity as Hell 4686, (b) a second absorption 
component which varies in strength and radial velocity like the C 1 1 1 
4647 emission, and (c) an emission core of constant strength whose 
radial velocity is equal to the system velocity (-25 km/sec). Hell 
4686 and Cl 1 1 4647 were chosen to define the radial velocity curves 
because they are the only strong, unblended lines which lie near the 
centers of the observed echelle orders. 

If we assume that the continuous spectrum of the primary is that 
observed at minimum light (justified by the disappearance of all 
spectral features varying in phase with the C 1 1 1 4647 emission at 
this time), and that all radiation from the primary which is incident 
upon the secondary is absorbed and re-emitted according to Planck's 
law, then we can use Tables I and III of Napier (1968) to reproduce 
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the first five light curves in Fig. 1, given the radii of the two 
stars relative to their separation, the inclination of the orbital 
plane to the plane of the sky, and the effective temperature of the 
primary. These tables take fully into account the effects due to the 
finite sizes and separations of the two components. The curves given 
in Fig. 1 were obtained using 0.12 and 0.18 for the radii of the 
primary and secondary, respectively, in terms of their separation; an 
inclination of 72° (for which the system just misses eclipse); and an 
effective temperature of 77,000 K for the primary. 

The effective temperature of 77,000 K for the primary was 
obtained using R, the ratio of the dereddened flux received by the 
short -wavelength IUE camera (1240 - 1945 A) to that received by the 
long-wavelength camera (1945 - 3120 A), at minimum light. 

Schonberner and Drilling (1984) have shown this quantity to be 
sensitive to the effective temperature, but to be relatively 
insensitive to interstellar reddening, surface gravity, and 
composition for n(He)/n(H) < 30%, and have calibrated R against 
effective temperature using both model atmospheres and stars whose 
effective temperatures have been determined by other methods. The 
goodness of fit of a 77,000 K model to the dereddened energy distri- 
bution of LSS 2018 near minimum light may be judged from Fig. 3 of 
Schonberner and Drilling (1984). The corresponding angular radius for 
the primary is l"5 x 10” 6 . 

Numerical experimentation shows that the amplitudes and shapes 
of the light curves are not very sensitive to the effective temper- 
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ature of the primary, but that they are extremely sensitive to the 
stellar radii, and that a change of only 5% in the radius of either 
star causes a noticeable departure from the observations. The light 
curves are not very sensitive to the inclination of the system for i 
> 53°, but cannot be fit at all for i <53°. The light curves alsn 
cannot be fit if Kurucz's (1979) modei atmospheres with log g = 4.5 
and solar composition are substitjtea for Planck's law in computing 
the spectrum of the secondary, as they predict a much larger Balmer 
jump than that which is observed. The use of these models is, 
however, unrealistic, as they assume that the atmosphere of the 
secondary is being illuminated from the inside rather than from the 
outside, as is the actual case. It would therefore be of interest to 
compute models with the appropriate boundary conditions, in order to 
try to predict the line spectrum and to see why the continuous 
spectrum of the secondary obeys Planck's law so well. 

The solid curve fit to the points representing the C 1 1 1 4647 
emission line strength in Fig. 1 was obtained by scaling the curve 
for the B magnitudes down to the observed amplitude for the Cl II 4647 
emission. The dotted curve is the best fit made on the assumption 
that the hemisphere of the secondary facing the primary is uniformly 
illuminated by the Cl 1 1 emission. It is seen that the second curve 
gives a better fit to the observations than the first. If we 
decrease the inclination of the system from 72° to 53°, the strength 
of the C 1 1 1 emission predicted by the dotted curve at minimum light 
increases to 20% of its strength at maximum light. Since neither 



- 10 - 


this nor an eclipse is observed, we conclude that the Inclination of 
the system lies between 53° and 72°. 

A complication arises in analyzing the radial velocity curve of 
the secondary because, according to the model employed, the center of 
light for the Cl 1 1 emission and the center of mass of the secondary 
do not coincide at quarter phase. However, it can easily be shown 
that the separation of these two points is 0.424 times the radius of 
the secondary if the hemisphere facing the primary is uniformly 
illuminated by the Cl 1 1 emission. If we assume synchronous rotation 
for the secondary, we then have from the amplitudes of the radial 
velocity curves: 

r b sin i = 361,000 km, and 
(r r - 0.424 R r ) sin 1 = 725,000 km, 
where r b and r r are the distances of the primary and secondary, 
respectively, from the center of mass of the system; R r is the radius 
of the secondary; and i is the inclination of the orbit. Application 
of Kepler's Third Law then gives 0.4 M Q (for the primary) and 0.2 M Q 
(for the secondary) as lower limits on the masses (i * 72°). The 
corresponding upper limits (for i * 53°) are 0.7 M 0 (for the primary) 
and 0.3 M 0 (for the secondary). The separation of the two components 
must lie between 0.008 and 0.010 AU, and the radii of the two stars 
are 0.2 - 0.3 R 0 (primary) and 0.3 - 0.4 R 0 (secondary). Comparison 
of the angular and linear radii of the primary yields 650 to 800 pc 
as the distance to the system, and an angular diameter of 180" for 
the planetary nebula (Drilling 1983) then gives 0.6 to 0.7 pc as its 
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linear diameter. The corresponding absolute magnitude for the 
primary (minimum light) is M y * 2.3 to 2.7. 

IV. DISCUSSION 

The reflection effect for LSS 2018 is one of the most pronounced 
such effects observed to date. The fact that the light variation can 
be predicted so well over such a large range in wavelength using the 
simplest possible assumption, that all radiation Incident on the 
secondary is absorbed and re-radiated according to Planck's law, is 
remarkable. It should be pointed out, however, that the masses 
derived above are, to a large extent, independent of this assumption. 
This is important, because LSS 2018 is the- only double-lined spectro- 
scopic binary known among the central stars, and because such strong 
constraints can be put on t.he inclination of this system, we have 
been able to make the first direct determination of the mass of a 
central star. 

If we plot LSS 2018 in Fig. 7 of Schonberner (1981), in which 
post-AGB evolutionary tracks for central stars are plotted in the M y - 
nebular radius plane (assuming a constant expansion velocity of 20 
km/sec), we find that the mass of the primary must be 0.57 M Q in 
order to be consistent with evolutionary theory. This is near the 
mid-point of the range in masses derived above, and our results thus 
support Schonberner' s conclusion that the mass distribution of 
centnl stars has a very sharp peak around C.58 solar masses. 

The masses and radii derived above for the secondary put it on 
or slightly above the mass-radius relation for main-sequence stars 
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given by Ritter (1982), which is based on data given by Popper 
(1980). As pointed out by Ritter (1982), this produces a problem In 
trying to explain the present state of the system as the result of 
common-envelope binary evolution, as the secondary has not had time 
to lose the required amount of mass and come back into thermal 
equilibrium since the ejection of the observed planetary nebula. 

According to the formula given by Kraft, Mathews, and Greensteln 
(1962), the period of LSS 2018 should be decreasing at the rate of 1- 
3 x 10" 5 days per century due to gravitational radiation. This will 
cause the secondary to fill its Roche lobe in 1-3 x 10 10 years. Since 
the primary will evolve into a white dwarf in a much shorter time 
than this, a cataclysmic binary will be produced. This may happen in 
a much shorter time if the system can find an additional means of 
losing angular momentum, such as the magnetic braking discussed by 
Verbunt and Zwaan (1981). Since the phase shift due to gravitational 
radiation will amount to 50-100 seconds in 500 years, the matter can 
eventually be settled by careful monitoring of the light variations. 
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FIGURE CAPTIONS 

Figure 1 . (a) 13.755 - V_. (b) 13.435 - B_. (c) 11.885 - _U. (d) 

2.5 log + 22.62, where is the total (dereddened) flux 
received between 1945 A and 3120 A. (e) 2.5 log f 2 + 21.21, 
where f 3 is the total (dereddened) flux received between 
1240 A and 1945 A. (f) 2.5 log (f 3 /f 4 + 1) + 0.25, where f 3 
is the total flux of Cl 1 1 4647 emission received and f 4 is 
the total flux received in a 3.4 A band of underlying 
continuum at minimum light. The squares represent observa- 
tions made from the ground and the triangles, IUE observa- 
tions. Open and closed symbols represent observations made 
on different nights or during different IUE shifts. Crosses 
represent very accurate photometric observations made using 
the B_ filter alone (averages over 1000 second intervals and 
two different nights). The curves are explained in the 
text. 

Figure 2 . (a) Heliocentric radial velocities for LSS 2018 given by 
the C 1 1 1 4647 emission line (squares) compared to a sine 
curve of amplitude 148 km/sec and phase shift n with respect 
to maximum light, (b) Heliocentric radial velocities for 
LSS 2018 given by the Hell 4686 absorption line (triangles) 
compared to a sine curve of 74 km/sec amplitude and zero 
phase shift. The system velocity is assumed to be -25 
km/sec in both cases. Open symbols represent observations 
made on April 12, and closed symbols, April 13. 
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